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We extend the program initiated in [T. Werlang et a/., Phys. Rev. Lett. 105, 095702 (2010)] 
in several directions. Firstly, we investigate how useful quantum correlations, such as entanglement 
and quantum discord, are in the detection of critical points of quantum phase transitions when the 
system is at finite temperatures. For that purpose we study several thermalized spin models in the 
thermodynamic limit, namely, the XXZ model, the XY model, and the Ising model, all of which 
with an external magnetic field. We compare the ability of quantum discord, entanglement, and 
some thermodynamic quantities to spotlight the quantum critical points for several different tem- 
peratures. Secondly, for some models we go beyond nearest-neighbors and also study the behavior 
of entanglement and quantum discord for second nearest-neighbors around the critical point at finite 
temperature. Finally, we furnish a more quantitative description of how good all these quantities 
are in spotlighting critical points of quantum phase transitions at finite T, bridging the gap between 
experimental data and those theoretical descriptions solely based on the unattainable absolute zero 
assumption. 
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I. INTRODUCTION 

Quantum phase transitions (QPTs) theoretically oc- 
cur at absolute zero temperature (T = 0) due to abrupt 
changes in the qualitative properties of the ground- 
state of a many-body system while varying its Hamil- 
tonian pQ. The tuning parameter (the quantity be- 
ing changed in the Hamiltonian) can be, for example, 
an external magnetic field or a coupling constant. By 
properly tuning the Hamiltonian one can reach a spe- 
cial point, the critical point (CP), where the ground 
state of the system undergoes a radical change, which 
strongly affects the macroscopic properties of the sys- 
tem. In the vicinity of a CP, a tiny change in the 
tuning parameter will favor one phase over the other. 
Since at T = there are no thermal fluctuations, these 
tiny changes are caused only by quantum fluctuations, 
whose origin can be traced back to the Heisenberg uncer- 
tainty principle. Some well-known examples of QPTs are 
the paramagnetic-ferromagnetic transition in some met- 
als [2], the superconductor-insulator transition [3 , and 
superfluid-Mott insulator transition [4]. 

Although reaching the absolute zero is impossible since 
it would violate the third law of thermodynamics, QPTs 
can be detected if one reaches temperatures near the 
absolute zero. More precisely, when the system's de 
Broglie wavelength is greater than the correlation length 
of the thermal fluctuations one achieves a temperature 
low enough to see a QPT. In this regime thermal fluc- 
tuations are unable to excite the system from its ground 
state and a phase transition solely driven by quantum 
fluctuations is still possible to be seen. 

Many theoretical tools employed in the calculation of 
these CPs are usually subjected to the T = restric- 
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tion. Interesting examples are the behavior of bipartite 
[5], multipartite [6], and generalized entanglement [7] for 
spin chains at T = 0. Quantum discord, a measure of 
all quantum correlations present in a system [8, 9 , has 
also been widely used to characterize the CPs at T = 
[T0l [TT] . The extremal values and the behavior of the 
derivatives of these quantities can signal the existence 
of a CP without the knowledge of an order parameter, 
i.e., without the knowledge of what macroscopic quantity 
abruptly changes during the QPT. 

All experiments involving QPTs are, nevertheless, per- 
formed at low temperatures, not at absolute zero, which 
forbid a direct connection between the measured data 
and those theoretical results developed for T = 0. Re- 
cently, however, we have studied, in the thermodynamic 
limit, the XXZ Hamiltonian with no magnetic field and 
shown that quantum discord (QD) is able to detect CPs 
at finite T while other thermodynamic quantities and en- 
tanglement fail in this task [12]. For that model we have 
shown that the thermal quantum discord (TQD), i.e., 
QD computed for a system described by the canonical 
ensemble, possessed a characteristic behavior at the CP 
robust enough to survive an increasing T and, therefore, 
spotlight a CP of a QPT at finite T. 

Our goal in this work is to extend the study initiated in 
[T2] in at least two directions. We want, on one hand, to 
study different models to see whether or not TQD is still 
a good CP detector at finite T. We want to answer the 
following question: Does TQD possess a clear distinctive 
behavior at the CP strong enough to be seen at T > 0? 
On the other hand, we want to be more quantitative 
in determining the ability of TQD, and other quantities 
such as entanglement, to correctly detect the CP of a 
QPT at finite T. How close to the actual CP are the 
estimated CPs obtained at finite T via TQD? How do 
these estimated CPs at T > approach the correct CP 
as we decrease T? Those are the questions we also want 
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to answer in this paper. 

To achieve these goals, we consider here the XXZ 
model in the presence of an external magnetic field and 
the XY model in a transverse field. For both models we 
analyze the behavior of TQD and entanglement for T > 
around their CPs. For the XXZ model with an external 
field we compare the efficiency of TQD with respect to 
entanglement in the detection of CPs as the temperature 
increases and also with other thermodynamic quantities. 
For the XY model we study the behavior of TQD and the 
entanglement of its first and second nearest-neighbors as 
a CP detector. 

In order to make this work self-contained, we organize 
its presentation as follows. In Sec. II we review some 
concepts from information theory necessary to the con- 
struction of two measures of quantum correlations: QD 
and entanglement of formation (EoF). In Sec. Ill we dis- 
cuss qualitatively and quantitatively the role of quantum 
correlations as a tool to detect QPTs at zero and finite 
temperatures in the XXZ and XY models, both with ex- 
ternal magnetic fields. Our final remarks and conclusions 
are presented in Sec IV. 

II. QUANTUM CORRELATIONS 
A. Quantum Discord 

In the context of classical information theory [13] , cor- 
relation is a measure of the dependence between two or 
more random variables or observed data values. The to- 
tal correlation between two random variables A and B 
with probability distribution p a and p b , respectively, is 
given by the mutual information (MI) [13] 

X 1 (A : B) = H(A)+H(B) - H(A,B), (1) 

where H(X) — — ^2 k p x ^og 2 p x is the Shannon entropy 
with p x the probability distribution of the random vari- 
able X. If p a , 6 denote the joint probability distribution 
of the variables A and B then the joint entropy H(A,B) 
can be written through the Bayes' rule p a \ b = p a ,b/Pb [2] 
as 1-L(A,B) = H(A\B) +U{B), where the conditional 
entropy H(A\B) = - E a ,6^ lo S2Pa|6 quantifies how 
much uncertainty is left on average about A when one 
knows the values of B. Note that the conditional entropy 
is always non-negative and not necessarily a symmetric 
quantity, i.e., H(A\B) may be different from %{B\A). 
Using this relation between H(A, B) and H(A\B) we can 
obtain an equivalent version for the MI as given in 0, 

1 2 (A:B)=H(A)-H(A\B). (2) 

For a quantum bipartite system AB described by the 
density operator p ABl a measure of their total corre- 
lation is obtained directly from Eq. ([I]) replacing the 
Shannon entropy by the von- Neumann entropy, S(X) = 
S(p x ) = -Tv(p x \og 2 p x ), i.e., X\(A : B) = S{A) + 
S(B) — S(A, B). However, a quantum analog of the sec- 
ond version of MI, Eq. (J2|, cannot be obtained so easily 



because Bayes' rule are not always valid in the quan- 
tum case [15 . For example, assuming the validity of this 
rule, the quantum conditional entropy would be defined 
as S(A\B) = S(A,B) - S(B). But this definition im- 
plies a negative value for S(A\B) if one computes it for 
the quantum state = ( 1 00) + |11)) /V%, showing that 
this quantity cannot be interpreted as in the classical 
case, where it is always positive. 

To ensure that the quantum conditional entropy has 
the same meaning as in the classical case, Henderson and 
Vedral [9] defined the conditional entropy as S q (A\B) = 
min{ Mb } ^ b p b S(A\B = b) such that the minimization 
is given over generalized measurements {M b } [13] . with 
^2 b M b = 1^, 1b the identity operator that acts on 5, 
M b > for all 6, and S(A\B = b) = S(p A \ b ) where 

pA\b = — (1a ® M b )p AB (l A <8> M 6 ), 
Pb 

with p b = Tr [(1 A <g> M b )p AB (l A <g> M b )}. The expres- 
sions S q (A\B) and S(A\B) coincide with H(B\A) for 
classical systems. Using S q (A\B) the second quantum 
version for MI, Eq. can be written as X 2 (A : B) = 
S(A) — S q (A\B). The difference between these two ver- 
sions of quantum mutual information, 

D{A\B)=J^{A:B)-ll{A:B), 

= S(B)-S(A,B)+S q (A\B), (3) 

was called quantum discord by Ollivier and Zurek [8 
and interpreted as a measure of total quantum corre- 
lation. The notation D(A\B) indicates that the mea- 
surements were performed on the subsystem B. When 
S(A) = S(B) we have D(A\B) = D(B\A) = D (p AB ) 
and therefore the QD is a symmetric quantity. As demon- 
strated recently [16] the state p AB has D(A\B) = if, 
and only if, there exists a von Neumann measurement 

{nf } such that PAB = Zj (u (8) nf ) PAB (i A nf ) 

[T7] , Therefore, a state p AB with D = has necessar- 
ily the form p AB = J2j,kPj,k \ ® \<t>k) (<f>k\ with 
^Z j}k Pj,k = 1 and {\ipj)} and {\4>k)} sets of orthogonal 
states. The states with D = are completely classically 
correlated in the sense of [18] . An example of such state is 
p AB = (|0) (0| |0) (0| + |1> (1| |1) (1|) /2. This result 
shows that the existence of quantum correlations is due 
to the principle of superposition that allows describing a 
system through a set of non-orthogonal states. Here, and 
analogously to entanglement computed for thermal states 
[19], we use the nomenclature coined in [20 , namely, 
thermal quantum discord (TQD), to refer to QD com- 
puted for states described by the canonical ensemble. 

The main difficulty in the determination of a closed ex- 
pression for QD lies in the complicated minimization pro- 
cedure for calculating the conditional entropy S q (A\B). 
For two-qubit systems the minimization over generalized 
measurements is equivalent to a minimization over pro- 
jective measurements (von Neumann measurements) [21]. 
In this case, the minimization procedure can be done nu- 
merically for general two-qubit states [22 while an ana- 
lytical expression can be obtained only for a subclass of 
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the so-called X-states (see Appendix). In this paper we 
employed both strategies, i.e., numerical minimization 
and, when available, closed expressions for computing 
QD. 



B. Entanglement 

Next, we introduce an important kind of quantum cor- 
related states, the entangled states. A quantum bipar- 
tite state described by pab is said to be entangled if, 
and only if, it cannot be written as a separable state 
Pab = ^ZjPjpf ® pf, where pj > and J2jPj = 1 
[23| [24] . Although there are states with non-zero quan- 
tum discord created only by local operations and clas- 
sical communication (LOCC), entanglement cannot be 
generated in this way [25] . For pure states all quantum 
correlated states are entangled. The situation is more 
complicated in the case of mixed states. For example, 
the Werner state p w = a \ + (1 — a) 1/4, with 

|^-) = (|01) - 1 10>) /y/2 and a e [0, 1], is separable for 
a < 1/3 and non-local (violates a Bell-like inequality) for 
a > l/y/2. However, the Werner state have null quantum 
discord only for a = 0. 

In this paper the measure of entanglement used is the 
Entanglement of Formation (EoF) [25 , which is defined 
as EoF(pab) = min^. PjS(pA(B)), where the mini- 
mization is over all ensembles of pure states {pj, \ipj)} 
such that pab = ^ZjPj l^j) (V^'l an d S(pa(b)) is the von 
Neumann entropy of the reduced state of either A or 
Pa(b) — T^b(A) (pab)- The EoF quantifies, at least 
for pure states and asymptotically, how many singlets are 
needed per copy to prepare many copies of pab using only 
LOCC. It is worth emphasizing that it coincides with the 
quantum discord when one deals with pure states. An an- 
alytic closed expression for EoF is given in the Appendix. 



III. RESULTS AND DISCUSSIONS 
A. The XXZ Model 

The first model tackled in this paper is the 
one-dimensional anisotropic spin- 1/2 Heisenberg chain 
(XXZ) subjected to an external magnetic field in the z- 
direction. The Hamiltonian of this model can be written 



h xxz = jJ2 + + 

u L 



(4) 



where a? (a = x, y, z) are the usual Pauli matrices acting 
on the j-th site and crg +1 = erf (periodic boundary condi- 
tions). Here h is the external magnetic field, J is the ex- 



change constant ( J = 1) and A is the anisotropy param- 
eter. We are going to consider only the case A > 0, how- 
ever the negative values A < can be obtained by revers- 
ing the sign of J followed by a canonical transformation 
[27] , The critical regime |A| < 1 can be parametrized by 
A = cos 7 and the non-critical regime |A| > 1 is given 
by A = cosh 7. The thermodynamical properties are 
obtained through the free-energy / = — 
where Z = Tr {exp (— j3H xxz )} is the partition function, 
with j3 = 1/kT and the Boltzmann constant k is set to 
unity. 

The free-energy is calculated by Bethe ansatz tech- 
niques and can be written as follow, 



f = e -^{V*\nBB) (0), 
where the ground state energy eo is given by 

) S in7 joo sinh(^-l^ fc 



(5) 



e J cos 7 
7 



00 2 sinh : 



cosh 7 — 2 sinh 7 



- cosh \ 
-00 cosh 7/c ' 



dk, < A < 1, 



A > 1, 



and 



V(x) = jgS,™ 



< A < 1, 
A > 1. 



(6) 



(7) 



jk= — oo 2cosh7fc' 

The symbol * denotes convolution / * g(x) — J__ a f(x — 
y)g(y)dy, where a— ^ooforO<A<l and a = tt/2 for 
A > 1. 

The auxiliary functions b(x), b(x), and its simply re- 
lated functions B{x) — b(x) + 1 and B(x) = b(x) + 1 are 
solution to the following set of non-linear integral equa- 
tions [28] . 

In 6(a) = d+(x) + (K*lnB)(x) - (if* In B) (x + i 7 ),(8) 
In 6(a) = d-(x) + (K*]nB)(x) - (K*lnB)(x - i 7 ).(9) 
The driving term d±(x) is given by 



Sin 7 7T 



d±(x) 



-2 J/3 



-2J/3sinh 7 ^ 



/3h 7T 



7 cosh (7^/7) 2 7T— 7 ' 



< A < 1, 



k— — 00 cosh 7 A; 



±!f, A>1, 



and the Kernel function 



(10) 



roo sinh (-7T — 27)-? 



K(x) 



Eoo 
k= 



2 sinh (7T— 7) cosh 
7|fc| A2kx 



T dk, 



-co cosh 7/c 



< A < 1, 
A > 1. 



(11) 

The system density operator p is described by the 
canonical ensemble, i.e., p = exp (— f3H xxz )/Z. The 
nearest-neighbor two-spin state is obtained by tracing all 
but the first two spins, pi 5 2 = Tr£_2(p). The Hamilto- 
nian Q exhibits both translational invariance and U(l) 

) , therefore the reduced 



invariance ^ H xxz ,Y,j 
density matrix will be given by 


1 P22 P23 
4 



Pl,2 



/ Pu 






P23 P22 






(12) 
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where 

Pll = 1 + 2 <<7*) + o\) , 
P22 = 1 - (0"Jcr|) , 
P44 = 1 " 2 (O + (<7f<7|) , 
P23 = 2 Kerf) . 

The magnetization and the two-point correlations above 
are obtained in terms of the derivatives of the free-energy 

EH 



(a z ) = -2d h f, 
(a]a] +1 )=d A f/J, 



u - Ad A f + h (a z ) 
2J : 



(13) 
(14) 

(15) 



<*W + i> = <°W + i> = H± ^. A = 1 > (16) 

where u = dpjPf) is the internal energy. 

Using ( fl2| ) we can obtain analytical expressions for 
TQD and EoF between two nearest-neighbor spins 
through Eqs. ( |A.3| ) and (A.4), respectively. See Ap- 
pendix for details. 

As discussed in [30 j, the XXZ model has two CPs in 
the absence of external field (h = 0). An infinite-order 
transition at A = 1, with the ground state changing from 
an XY-like phase ( — 1 < A < 1) to an Ising-like antifer- 
romagnetic phase for A > 1, and a first-order transition 
at A = —1, from a ferromagnetic phase (A < —1) to 
the critical antiferromagnetic phase (— 1 < A < 1). In 
[lOj [TT] it was shown that both quantum discord and en- 
tanglement are able to detect these CPs associated to 
QPTs at T = 0. In order to extend these zero tempera- 
ture studies of the XXZ model, we take into account the 
action of the external field h. The effect of the external 
field is to shift the critical points to higher values of A 
[31]. The critical point associated to the infinite-order 
transition A in f is determined by the following equation 



h = 4Jsinh(^) 



oo 

E 

I— — oo 



(-1)" ; 

cosh(n^) ' 



(17) 



with rj = cosh 1 (Aj„/), while the critical point associ- 
ated to the first-order transition Ai by the equation 



/i = 4J(l + Ai). 



(18) 



Table I shows the critical points associated to first- and 
infinite-order QPTs for different values of the external 
field h. 

We begin our analysis by studying the behavior of 
TQD and EoF as a function of the tuning parameter 
A for T = and h > 0. These results are shown in Fig. 
[I] For h = |T0j [11] the infinite-order transition is char- 
acterized by a global maximum of TQD and EoF at the 
CP Ai n f — 1 ( see Pig- [l^)- But, besides being a max- 
imum, TQD also has a cusp in this CP. This behavior 
implies that the first derivative of TQD is discontinuous 
at Ai n f = 1 while the second derivative is divergent. For 



TABLE I. Critical points associated to the first-order QPT 
(Ai) and to the infinite-order QPT (Ainf) for different values 
of the external field h. 



h = 



h- 



h = 6 h=12 



Ai (first-order) -1.00 -0.50 0.50 2.00 

A^f (infinite-order) 1.00 2.15 3.30 4.88 



the first-order CP Ai = — 1 , it was shown in [TT] that 
for h = the EoF and the TQD are discontinuous in the 
CP. 

However, for h > Figs. [I}>d show that the infinite- 
order transition is characterized through TQD and EoF 
by a cusp at the CP rather than a global maximum or 
minimum. For a different spin model it has been recently 
shown [32] that the extreme points of entanglement does 
not necessarily indicate an infinite-order transition. Here 
we show an example where an infinite-order CP is pointed 
out not by a global maximum/minimum of the value of a 
quantum correlation (EoF or TQD). Rather, it is a dis- 
continuity in its first derivative that spotlight the infinity- 
order QPT. Furthermore, the first-order CP Ai is well 
determined by TQD and EoF, as can be seen for h = 6 
(Fig. [l^) and h — Yl (Fig. [l]i). Both quantities are zero 
for A < Ai and nonzero for A > A x while their first 
derivatives are divergent at Ai. It is worth mentioning 
that the cusp-like behavior in Fig. [TJbetween the two CPs 
is due to an exchange in the set of projectors that min- 
imizes the quantum conditional entropy [T2] , Eq. (A.2), 



and so far could not be associated with any known QPT 
for this model. 
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FIG. 1. (color online) QD (black solid line) and EoF (red 
dashed line) as functions of the tuning parameter A and ex- 
ternal field h for the XXZ model in the thermodynamic limit 
at T = 0. (a) h = 0,(b)h = 2, (c) h = 6 : and (d) h 12. The 
dashed black vertical lines denote the CPs Ai and A in f. Here 
and in the following graphics all quantities are dimensionless. 



The next step is to examine how the characterization 
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of CPs through TQD and EoF is affected by tempera- 
ture when ft > 0. To this end, let us start looking at 
the behavior of TQD and EoF for finite T and ft = 0. 
In Figs. [2^i and (2}d we plot, respectively, TQD and EoF 
as a function of A for some values of T > and ft = 0. 
As discussed in our previous work [T2] . and illustrated 
in Fig. [5J for ft = the first-order derivative of TQD 
is discontinuous at the CP A = 1 not only at T = 0, 
but also at T > 0. On the other hand, the EoF is max- 
imum at this CP only when T = and as T increases 
the maximum moves to higher A values (A > 1). The 
behavior of TQD in this case is related to the minimiza- 
tion process of the quantum conditional entropy (A.2). 



We observe that for finite T there is an exchange in the 
functions that minimize the conditional entropy which in 
turn leads to a sudden change in TQD [33], character- 
ized by a discontinuity in its first-order derivative. It is 
worth noting that for ft = we found a simple closed 
expression for TQD involving the two-point correlation 
functions. This closed expression together with the fact 
that at A = 1 we have = |(cr|cr| +1 )| allowed 

us to explain the sudden change occurring at this QPT 
forT>0[l2]. 




FIG. 2. (color online) (a) TQD and (b) EoF as a func- 
tion of A for ft = and, from top to bottom when A < 1, 
kT = 0.02,0.1,0.5,1.0,2.0. The sudden change of TQD re- 
mains in the CP A^ n / = 1 as the temperature increases while 
the maximum of EoF is shifted to the right. 

Moving to the cases where ft > we keep the same 
set of parameters of Fig. [2] but the value of the exter- 
nal magnetic field. In Figs. [3] and [4] we use ft = 6 and 
ft = 12, respectively. Looking at Figs. [3] and [4] we note 
that as the temperature increases the cusp-like behavior 
of TQD is smoothed while both the maximum of TQD 
and EoF decrease. Also, both curves of TQD and EoF 
become more smooth and broadened, turning these func- 
tions different iable in the CP. Interestingly, however, if 
the temperature is not too high the derivatives of these 




FIG. 3. (color online) (a) TQD and (b) EoF as a function 
of A for ft = 6 and, for A « 2, from top to bottom kT = 
0.02,0.1,0.5,1.0,2.0. 




FIG. 4. (color online) (a) TQD and (b) EoF as a function 
of A for ft = 12 and, for A « 3, from top to bottom kT — 
0.02,0.1,0.5,1.0,2.0. 



quantities will still keep some features of their behavior 
at T = 0. Indeed, as we have shown for T = 0, the CP 
Ai is characterized by a divergence in the first derivative 
of both TQD and EoF while the CP A^ n f by a divergence 
in the second derivative. As the temperature increases, 
though, these divergences disappear but the appropriate 
derivatives of TQD and EoF assume their maximum val- 
ues around the CPs, a property resembling the behavior 
at T = 0. To illustrate this point we plot in Fig. [5] the 
derivatives of TQD for ft = 12 as a function of A and 
for kT = 0.02,0.1,0.5. Note that the first-order deriva- 
tive of TQD has a local maximum near the CP Ai = 2 
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FIG. 5. First-order derivative of TQD as a function of A for 
(a) kT = 0.02, (b) kT = 0.1, and (c) kT = 0.5. Second-order 
derivative of TQD as a function of A for (d) kT = 0.02, (e) 
kT = 0.1 and (f) kT = 0.5. The maximum of the first and 
second derivatives of TQD are very close to the CPs Ai = 2 
and Ai n f = 4.88, respectively. Here we fixed ft — 12. Dashed 
vertical bars indicate the CPs. 



for kT = 0.02,0.1,0.5. The infinite-order QPT is also 
well characterized by a local maximum in the second- 
order derivative of TQD near the CP A in f — 4.88 for 
kT = 0.02,0.1,0.5. This procedure to determine, or at 
least estimate with great accuracy, the CPs with finite 
temperature theoretically computed data has been tested 
for other values of ft and it was proved valid. The same 
analysis can be applied to EoF but, as we show in the 
following, TQD gives better results than EoF for all sets 
of parameters. 

In addition, in Fig. [6] we show, for several temper- 
atures, many thermodynamic quantities for the infinite 
spin chain and also the pairwise correlations as a func- 
tion of A for ft = 12. For low temperatures the CPs can 
be estimated by some thermodynamic quantities such as 
the magnetic susceptibility, entropy, and specific heat. 
Although TQD is not the only quantity that can deter- 
mine or estimate a CP close to T = (for ft > 0) it is the 
best CP detector compared with other quantities tested 
here. This is shown in Fig. [7| where we plotted the CPs 
determined via TQD, EoF, and the pairwise correlations 
(erf af) and (erf af) for several temperatures. We chose 
to work with ft = 6 and ft = 12 but similar results are 
obtained for other values of the magnetic field. 

The procedure adopted to estimate the CPs when T 
is finite comes from the behavior of TQD, EoF, and the 
two point correlations near or at the CPs. For example, 
in Figs. [2j|3j and [6] we see that in the vicinity of the CPs 
these quantities, as a function of A, are either maximum 
or possess discontinuous/divergent first /second deriva- 
tives. Therefore, by computing the appropriate deriva- 
tive and looking for its maximum/minimum we can esti- 
mate the value of the CPs even for finite T. Obviously, 
for T = these extrema are exactly at the CP. As can 




FIG. 6. (color online) Thermodynamic quantities for the 
XXZ model in the thermodynamic limit for ft — 12 and kT = 
0.02 (black line), kT = 0.1 (red line), kT = 0.5 (blue line), 
kT =1.0 (green line), and kT = 2.0 (orange line). Here the 
CPs are Ai = 2 and A inf = 4.88. For A < 1 and for the 
curves of the magnetization and the two-point correlations 
kT increases from top to bottom while the opposite happens 
to the curves of the other quantities. 



be seen in Fig. [7| the CPs estimated by TQD are closer 
to the values given in Table I than those obtained using 
EoF or the pairwise correlations. Indeed, from zero to 
kT ~ 1.0 the extrema of TQD is closer to the actual CP 
than those from the other quantities (see Fig. [7]). 

To complement our studies about the XXZ model we 
present in Fig. [8] the behavior of TQD and EoF as a 
function of temperature for ft = 0. First, we observe 
that EoF tends to disappear suddenly as the tempera- 
ture increases while TQD is always nonzero. This sudden 
disappearance may be retarded by increasing the value 
of the anisotropy A. The insets in Fig. [8] show a "re- 
growth" [20 of TQD and EoF with temperature. This 
interesting feature was already seen for the behavior of 
TQD of a two-spin XXZ model with no field [20 although 
the regrowth does not show up for EoF [34]. Therefore, 
the results here presented show that the regrowth of EoF 
with temperature, while not observed for two-spin chains 
without the presence of an external field [34 , is possible 
in the thermodynamic limit. 



B. The XY Model 



In this section we discuss the ability of TQD and EoF 
to point out CPs of a QPT in the infinite one-dimensional 
XY model in a transverse field for finite T. The Hamil- 
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tonian that describes this model is given by 



0.0 0.2 0.4 0.6 0. 

kT 



1.0 0.0 0.2 0.4 0.6 0.8 1.0 

kT 



FIG. 7. (color online) CPs determined by TQD (square), 
EoF (circle), (erf erf) (up arrow), and (afa^) (down arrow) 
as a function of T. In (a) and (b) we have h = 6 with a 
first (a) and an infinite-order (b) CP; in (c) and (d) we have 
h = 12 with a first (c) and an infinite-order (d) CP. Dashed 
horizontal lines indicate the CPs. To determine the CPs at 
finite T we used the derivatives of TQD and EoF. As stated 
in the text, for some temperature ranges, the derivatives of 
these quantities at finite T retain some information about the 
divergences or discontinuities that exist only at T = 0. Note 
that the estimated CPs given by TQD is always closer to the 
exact value than the CPs coming from other quantities. 




FIG. 8. (color online) TQD (bottom) and EoF (top) as a 
function of T for different values of the tuning parameter A. 
The insets show the regrowth of TQD and EoF. For A < 1 
(solid curves) we note that as we increase A we also increase 
TQD and EoF. For A > 1 (dashed curves) this only happens 
for EoF at high temperatures. 



H 



A 



xy 



2 E[( 1 +«i + ( 1 -«i] 



L 



(19) 



The strength of the inverse of the external transverse 
magnetic field is represented by A while 7 provides the de- 
gree of anisotropy. For 7 = ±1 we obtain the transverse 
Ising model while 7 = corresponds to the XX model in 
a transverse field [35 j, which is the A = XXZ model 
discussed before with J < 0. As is well-known the XY 
model undergoes a second-order QPT (Ising transition 
[36] ) at the CP A c = 1 [10] , that separates a ferromag- 
netic ordered phase from a quantum paramagnetic phase. 
For A > 1 one can observe another second-order QPT 
at the CP 7c = called anisotropy transition [35] [37] . 
Differently from the Ising transition, that is due to the 
action of the external field, this transition is driven by 
the anisotropy parameter 7 and separates a ferromagnet 
ordered along the x direction and a ferromagnet ordered 
along the y direction. Although both QPT have the same 
order they belong to different universality classes [35] [37] . 

The XY Hamiltonian is ^-symmetric and can be ex- 
actly diagonalized [35] in the thermodynamic limit L — >• 
00. Due to the translational invariance the two-spin den- 
sity operator pij for spins i and j depends only on the 
distance k = \j — i\ between them and thus pij = po,fc- 
Therefore, the reduced density operator po,/c f° r the XY 
model at thermal equilibrium is [38] 



P0,k = 1 [Io,k + (cr Z ) (erg + <T%)] 

1 ^ \ ^ /^Q' ck\ ^ol ex 
a=x,y,z 



(20) 



where io,fc is the identity operator of dimension four. The 
transverse magnetization (<j|) = (<j z ) is 



(1 + A cos <j)) tanh (flu^) 



d(j) 



(21) 



with uo^ = (7A sin <j)) 2 + (1 + A cos (j)) 2 /2. The two- 
point correlation functions are given by 



/ ^.x ^x\ 



G-i G-2 
Gq G-i 



Gk- 


-2 Gk-3 




G_! 


Gi 


Go 


G_ 


fe+2 


G 2 


Gi •• 


G_ 


fe+3 



Gfe Gfe_i 



G-k 
G-k+i 



Gi 



(22) 

(23) 
(24) 
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where 
Gk 



tanh (PoJcf)) cos (k(j))(l + A cos (j>) 



7 A 



/ tanh (Puq) sin (kef)) sin 
Jo 



def) 



(Ob and (TO) of the 



Using these relations and Eqs. 
Appendix, EoF and TQD can be computed analytically 
for this model. 

The ability of TQD and EoF to detect a CP for this 
model at T = was studied in [39 , where TQD and 
EoF from first to fourth nearest-neighbors was computed. 
They showed that TQD between far neighbors can still 
characterize a QPT while EoF is zero. This result shows 
an advantage of the quantum discord to the entangle- 
ment in the detection of CP. See also [40] for the behavior 
of entanglement and TQD for small T in the symmetry 
breaking process of the XY model. Our goal here is to 
determine whether the quantum discord for the first and 
second neighbors has some advantage with respect to en- 
tanglement to spotlight a CP at finite T. The procedure 
adopted here is similar to the one employed in the anal- 
ysis of the XXZ model, i.e., we compute the appropriate 
derivatives around the CP and look for its extremum val- 
ues as indicators of QPTs. 

In Figs. |9| and [l0| we plot, respectively, TQD and EoF 
for first nearest-neighbors and second nearest-neighbors 
as a function of A for kT = 0.01,0.1,0.5 and 7 = 
0,0.5,1.0. It is important to note that TQD is more 



(a) 









o 0.1 




FIG. 9. (color online) (a)-(c) TQD and (d)-(f) EoF as a 
function of A for kT = 0.01 (black/solid line), kT = 0.1 
(red/dashed line) and kT = 0.5 (blue/dotted line) for nearest- 
neighbors. We use three values of 7 as shown in the graphs. 

resistant to thermal effects than EoF. For example, for 
kT = 0.5 EoF for the first-neighbors is zero or close to 
zero for almost all A while TQD is always non-null. For 
second-neighbors the situation is more drastic, EoF is al- 
ways zero for kT = 0.5. The previous discussion also 



(a) 


y=o ■ 








(b) 


y=o.5 




(c) 


Y=1 








FIG. 10. (color online) (a)-(c) TQD and (d)-(f) EoF as a 
function of A for kT = 0.01 (black/solid line), kT = 0.1 
(red/dashed line) and kT = 0.5 (blue/dotted line) for sec- 
ond nearest-neighbors. We use three values of 7 as shown in 
the graphs. 



highlights a qualitative difference between the behavior 
of quantum correlations as measured by quantum discord 
or entanglement, especially for 7 = 0.5. 

Since the Ising transition is a second-order QPT, the 
critical point at T = is characterized by a divergence 
or discontinuity in the first derivative of TQD or EoF 
[T0l fTTj ; and if the first derivative is discontinuous then 
the CP can also be found through the divergence of the 
second derivative. Similarly to the XXZ model studied 
in Sec. Ill A[ the singular behavior of TQD and EoF at 
the CP is attenuated as the temperature increases. How- 
ever, one can still extract useful information concerning 
the CPs until a certain temperature. Beyond this tem- 
perature the behavior of the derivatives of TQD or EoF is 
not able to determine the CPs without ambiguity. With 
that in mind, we use the same strategy of Sec. |III A to 
determine/estimate the CPs with theoretically computed 
data at finite T: if the first derivative of TQD or EoF is 
divergent at T = then the CP is pointed out by a local 
maximum or minimum at T > 0; if the first derivative 
is discontinuous at T = then we look after local max- 
imum or minimum in the second derivative for T > 0. 
The results of our analysis are shown in Fig. 



11 where 



the CPs estimated by TQD and EoF for first and second 
neighbors are plotted as a function of kT for 7 = 0,0.5, 
and 1.0. 

Looking at Fig. 11 we note that the CP of the Ising 
transition is best estimated by TQD than EoF. While 
for 7 = 0.5 and 7 = 1 the first- neighbor TQD is only 
slightly better than EoF to estimate the CP , for 7 = 
the CP computed by TQD is clearly closer to A c = 1 than 
the one estimated from EoF. Moreover, the temperature 
below which TQD is a good CP detector is higher for 
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FIG. 11. (color online) CPs determined by TQD (black- 
squares) and EoF (red-circles) as a function of kT for (a)- 
(c) first nearest-neighbors and for (d)-(f) second nearest- 
neighbors. Dashed horizontal lines indicate the CPs. The 
values of 7 used were 0, 0.5 and 1.0, as shown in the graphs. 
The derivatives of TQD and EoF were used to estimate the 
CPs as explained in the text. Note that at panel (b) the 
curves for TQD and EoF coincide and that at panels (d) and 
(f) we only have two points for EoF since for greater T it is 
zero. 



7 = than for 7 = 0.5 or 7 = 1 (see Fig. [TT^i) . When 
we look at the second- neighbors TQD and EoF, we see 
that TQD is still better than EoF to estimate the CP. 
And in particular for 7 = 1.0, EoF is only useful for very 
tiny T while TQD works relatively well for higher values 
of T. It is worth noting that for small T TQD seems to 
approach the actual CP almost linearly as we decrease T. 
Therefore, a relatively small number of points for a couple 
of small temperatures allows one to correctly extrapolate 
to the exact value of the CP at T = 0. This fact is also 
observed for the behavior of TQD in the XXZ model. An 
interesting future investigation, lying beyond the scope 
of this manuscript, would be to rigorously determine the 
functional form by which TQD and EoF approach the 
exact CP for decreasing T. 

We close this section studying the behavior of TQD 
and EoF near the critical point 7 C = 0, which is asso- 
ciated to the anisotropy transition that occurs only for 
A > 1. Hence, for definiteness, in what follows we set 
1.5. In Fig. [12] we plot TQD and EoF for the first 



A 

and second-neighbors as a function of 7 for five different 
values of temperature, kT = 0.001, 0.1, 0.5, 1.0, and 2.0. 



Looking at Fig. 12 we see that TQD and EoF are max- 
imum at the CP 7c = for first- and second-neighbors, 
keeping this feature as T increases. However, only TQD 
possesses a cusp-like behavior at the CP. Furthermore, 




-0.4 -0.2 0.0 0.2 0.4 -0.4 -0.2 0.0 0.2 0.4 



FIG. 12. (color online) Top panels TQD and bottom pan- 
els EoF as a function of 7. From top to bottom kT — 
0.001,0.1,0.5,1.0, and 2.0. Left panels are first nearest- 
neighbors and right panels second nearest-neighbors. Here 
we fixed A = 1.5 and the CP is j c = (indicated by dashed 
vertical lines). 



for first-neighbors TQD presents the same pattern (max- 
imum with a cusp-like behavior) up to kT = 2.0 while 
EoF vanishes for kT > 1.0. This situation is more dras- 
tic for second-neighbors: while TQD can detect the CP 
even for values near kT = 1.0, EoF is nonzero only at low 
temperatures (kT < 0.1). Finally, as in the XXZ model 
with no external field, there is an exchange in the set of 
projectors that minimizes the conditional entropy (A. 2) 



at the CP. This feature is characterized by a discontinu- 
ity in the first-order derivative of TQD and explains the 
cusp-like behavior of TQD at j c = 0. 



IV. CONCLUSIONS 

In this paper we have studied the ability of quantum 
correlations to spotlight critical points of quantum phase 
transitions when the system under investigation is away 
from the absolute zero, the experimentally unattainable 
temperature where a quantum phase transition sets in. 
We have mainly focused our analysis in two distinct types 
of quantum correlations, namely, the entanglement and 
the quantum discord between pairs of spins. 

We first concentrated our efforts in computing those 
quantum correlations for an infinite chain described by 
the XXZ model with an external field in the z-direction 
and in equilibrium with a thermal reservoir at temper- 
ature T. By tracing out all but two nearest-neighbors 
we were able to calculate their entanglement and quan- 
tum discord for several T. For this model we also evalu- 
ated several thermodynamic quantities as well. Among 
all these quantities, we have shown that quantum discord 
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is the best critical point detector at finite T, furnishing 
rather precise values for the two quantum critical points 
for this model, a first order and an infinite order phase 
transition. 

We then moved to study other models, such as the infi- 
nite chain XY model and the Ising model, both of which 
with external transverse magnetic fields. This time we 
computed the entanglement and the quantum discord for 
first and second nearest-neighbors, assuming the system 
in equilibrium with a thermal reservoir at temperature 
T. For the second order phase transition usually referred 
to as the Ising transition, we observed that quantum dis- 
cord is better suited in estimating the correct value for 
the quantum critical point. For the other quantum phase 
transition, the so called anisotropy transition (also sec- 
ond order), quantum discord and entanglement give ex- 
cellent results for low temperatures. However, as the 
temperature increases quantum discord outperform en- 
tanglement as a quantum critical point estimator. 

Summarizing, the previous discussions showed that for 
all the cases here presented, quantum discord was the 
best critical point estimator when one deals with finite 
temperature systems. Note that we never needed to know 
the order parameter of the phase transitions to use these 
quantum correlations as critical point detectors. More- 
over, we have also showed that using the quantum discord 
(our best critical point estimator), and with theoretically 
computed data for several T, we can with reasonable ac- 
curacy infer the actual value of the critical point at T = 
by looking at the behavior of these quantities when T 
approaches zero. Therefore, we have given strong indi- 
cations that entanglement and in special quantum dis- 
cord are useful quantities in estimating quantum critical 
points when the available data are away from the abso- 
lute zero. This fact alone shows that these quantities are 
an important tool in the study of quantum phase tran- 
sitions in realistic experimental scenarios, where one is 
always working at finite T and has perhaps no clue of 
the order parameter of the phase transition taking place. 
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Appendix: Closed expressions for QD and EoF 

In Sec. |III[ after tracing out all spins but two in an 
infinite chain, the two-qubit reduced density matrix is of 



X-form, 



Pab 



/Pll P14 \ 

#22 P23 

P23 P22 

V p!4 #44 / 



(A.l) 



It is easy to see that for this density matrix S(A) = 
S(B) = - Yh=i Pi lo &2 with f3i = p n + p 2 2 and f) 2 = 
P44 + P22, implying t hat Q D is a symmetric quantity. For 
the density operator ( A.l ) this minimization can be done 
analytically, as previously demonstrated in [26 , resulting 
in the following expression, 



S q (A\B) = mm{f3 1 T(6 1 ) 
where F(0) ■ 



(A.2) 

^log 2 ^-^log 2 ^and 



1 

Pi 
1 



\P22 ~ Pll| 



\P22 ~ PAA\ 



^4(^14 - P23) 2 + {fill ~ P44) 2 , 
^4(^14 + P23) 2 + {fill ~ P44) 2 . 



Using (A.2) in Eq. ([3| and calcu latin g the joint entropy 
S(A, B) of the density operator (A.l), the quantum dis- 
cord is obtained as 



D (pab) = Pi log2 Pi + ^2 X i log 2 X 3 



3 = 1 



+S q {A\B), 



(A.3) 



Ai 



A 2 



Pll + P44 + V (Pll - P44) 2 + 4 |p 14 | 



where Xj are the eigenvalues of (A.l), 

1 " 
2 
1 
2 

A3 = P22 + |/023 I 
A4 = P22 — I #23 1 



Pll + P44 - A/ (Pll - P44) 2 + 4 |/0i4 1 



When P14P23 7^ Eq. (A.2) provides the correct 



value for the quantum conditional entropy. However, for 
P14P23 = the correct value of the conditional entropy 
is not always given by this expression [22j [41]. In the 
latter case the minimization depends on only one single 
parameter and should be numerically carried out. For all 
the cases studied here, however, the numerical calcula- 
tions agreed with the analytical expression of quantum 
discord. 



by 



Finally, for the density operator ( A.l ) the EoF is given 

a 

EoF(p AB ) = -g log 2 g-(l-g) log 2 (l - g), (A.4) 

where g = (1 + y/l - C 2 )/2 and C = 2 max {0, Ai, A 2 } 
is the concurrence with Ai = \pi±\ — y/p22p33 and A 2 = 

I #23 I - VP11P44- 
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